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Abstract 

The semi-classical Bloch-Boltzmann theory is at the heart of our understanding 
of conduction in solids, ranging from metals to semi-conductors. Physical systems 
that are beyond the range of applicability of this theory are thus of fundamental 
interest. This is the case of disordered systems which present quantum interferences 
in the diffusive regime, i.e. Anderson localization effects. This is also the case, for 
example, of systems that present magnetic or electric breakdown when submitted 
to an electromagnetic field. These exceptions, for which a full quantum transport 
theory must be developed have been intensively studied in the past and are now 
well known. 

It appears that in quasicrystals and related complex metallic alloys another type 
of breakdown of the semi-classical Bloch-Boltzmann theory operates. This type of 
quantum transport is related to the specific propagation mode of electrons in these 
systems. Indeed in quasicrystals and related complex phases the quantum diffu- 
sion law deviates from the standard ballistic law characteristic of perfect crystals in 
two possible ways. In a perfect quasicrystal the large time diffusion law is a power 
law instead of a ballistic one in perfect crystals. In a complex crystal the diffusion 
law is always ballistic at large time but it can deviate strongly from the ballistic 
law at sufficiently small times. We develop a theory of quantum transport that 
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applies to a normal ballistic law but also to these specific diffusion laws and we de- 
scribe the behavior of conductivity that results from these specific laws. As we show 
phenomenological models based on this theory describe correctly the experimental 
transport properties. Ab-initio calculations performed on approximants confirm also 
the validity of this anomalous quantum diffusion scheme. This provides us with the 
first ab-initio model of conductivity in approximants such as the a-AlMnSi phase. 

Although the present chapter focuses on electrons in quasicrystals and related 
complex metallic alloys the concept that are developed here can be useful for 
phonons in these systems. There is also a deep analogy between the type of quan- 
tum transport described here and the conduction properties of other systems where 
charge carriers are also slow, such as some heavy fermions or polar onic systems. 
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1 Introduction 



Immediately after the discovery by Shechtman et al. [1] of quasiperiodic in- 
termetallics one major question was raised about the physical properties of 
phases with this new type of order. In particular, one expected that the elec- 
tronic and thermal properties could be deeply affected. Indeed the description 
of electrons or phonons in periodic phases rests on the Bloch theorem which 
cannot be applied to a quasiperiodic structure. Within a decade a series of new 
quasiperiodic phases and approximant were discovered and intensively stud- 
ied. These investigations learned us that indeed the electrons and the phonons 
properties could be deeply affected by this new type of order. 

The first quasiperiodic alloys AIMn where metastable and contained many 
structural defects. As a consequence they had conduction properties similar 
to those of amorphous metals with resistivities in the range 100-500 fiQcm. In 
1986 the first stable icosahedral phase was discovered in AlLiCu. This phase 
was still defective and although its resistivity was higher (800 fiQcm) it was 
still comparable to that of amorphous metals. The real breakthrough came 
with the discovery of the stable AlCuFe icosahedral phase, having a high 
structural order. The resistivity of these very well ordered systems where very 
high, of the order of 10 000 /jficm, which gave a considerable interest in their 
conduction properties. Within a few years several important electronic char- 
acteristics of these phases were experimentally demonstrated. The density of 
states in AlCuFe was smaller than in Al, about one third of that of pure Al, 
but still largely metallic. The conductivity presented a set of characteristics 
that were either that of semi-conductors or that of normal metals. In par- 
ticular weak-localization effects were observed that are typical of amorphous 
metals. Yet the conductivity was increasing with the number of defects just as 
in semi-conductors. Optical measurements showed that the Drude peak, char- 
acteristic of normal metals, was absent. In 1993 another breakthrough was the 
discovery of AlPdRe which had resistivities in the range of 10 6 /zf2cm [2,3,4]. 
This system gave the possibility of studying a metal-insulator transition in 
a quasiperiodic phase. There are still many questions concerning electronic 
transport in AlPdRe phases. One difficulty concerns the homogeneity and the 
quality of samples which are crucial for transport properties but are difficult 
to determine exactly. 

Since the discovery of Shechtman et al. [1] our view of the role of quasiperiodic 
order has evolved. For electronic or phonon properties of most known alloys it 
appears that the medium range order, on one or a few nanometers, is the real 
length scale that determines properties. This observation has lead the scientific 
community to adopt a larger point of view and consider quasicrystals as an 
example of a larger class. This new class of Complex Metallic Alloys contains 
quasicrystals, approximants and alloys with large and complex unit cells with 
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possibly hundreds of atoms in the unit cell. 

In this chapter we shall concentrate on "the way electrons propagate" in a 
quasicrystal or in a complex metallic alloy The main objective is to show 
that the non standard conduction properties of some quasicrystals and related 
complex metallic alloys result from purely quantum effects and cannot be 
interpreted through the semi-classical theory of transport. This is of great 
importance since the semi-classical Bloch-Boltzmann theory is at the heart 
of our understanding of conduction in solids, ranging from metals to semi- 
conductors. This new type of quantum transport is related to the specific 
propagation mode of electrons in these systems. Indeed in quasicrystals and 
related complex phases the quantum diffusion law deviates from the standard 
ballistic law characteristic of perfect crystals in two possible ways. In a perfect 
quasicrystal the large time diffusion law is a power law instead of a ballistic one 
in perfect crystals. In a complex crystal the diffusion law is always ballistic at 
large time but it can deviate strongly from the ballistic law at sufficiently small 
times. It is this specific character that provides a basis for the interpretation 
of the strange conduction properties of AlCuFe, AlPdMn and probably also 
for those of AlPdRe. 

This chapter is organized as follows. Part 2 is the most technical part but it can 
be skipped by readers not interested by mathematical aspects. Part 3 presents 
a detailed physical interpretation of anomalous diffusion and low frequency 
conductivity laws in crystals and quasicrystals. Part 4 presents evidence of 
anomalous diffusion in experimental quasicrystalline and approximant phases. 

In part 2 we give some definitions and present the mathematical relations that 
exist between the low frequency conductivity, including the dc conductivity, 
and quantum diffusion. We consider also the relaxation time approximation 
(RTA) that allows to treat the role of disorder on quantum diffusion and 
conductivity. We demonstrate general formulas for quantum diffusion and low 
frequency conductivity (within the RTA) in periodic and quasiperiodic models 
of potential. 

In part 3 we focus on the physical interpretation and consequences of the 
formulas derived in part 2. On a general ground we discuss the limitations of 
the RTA and the possibility of a metal-insulator transition. We apply this to 
a general theory of low frequency conductivity and metal-insulator transition 
in crystals and quasicrystals. 

In part 4 we present briefly the experimental transport properties of phases 
such as AlMnSi, AlPdMn and AlCuFe or AlPdRe. These experimental trans- 
port properties indicate a conduction mode which is neither metallic nor semi- 
conducting. For the a-AlMnSi phase, recent ab-initio computations are pre- 
sented, which confirm the existence of an anomalous diffusion and allow for a 
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semi-quantitative ab-initio computation of conductivity. Concerning AlCuFe 
and related quasiperiodic phases, which cannot be addressed by band struc- 
ture calculations, we present a phenomenological model. This model based on 
anomalous quantum diffusion provides a coherent interpretation of the strange 
electronic transport of these systems. 

We conclude by a short summary. We discuss also briefly the link with other 
problems such as phonons in quasiperiodic systems or electrons in heavy 
Fermions systems or in polaronic systems. 



2 Quantum formalism for electronic transport 



In this section we give some definitions and recall general properties of the con- 
ductivity. We present the mathematical relations that exist between the low 
frequency conductivity, including the dc conductivity, and quantum diffusion. 
We consider also the relaxation time approximation (RTA) [5] that allows to 
treat the role of disorder on quantum diffusion and conductivity. We demon- 
strate general formulas for quantum diffusion and low frequency conductivity 
(within the RTA) in periodic and quasiperiodic models of potential. 



2.1 Impulse response and analytical properties of the conductivity 



Let us consider a system, at thermodynamical equilibrium, submitted to an 
impulse of electric field 

E(t) = E5(t) (1) 



where S(t) is the Dirac function. The resulting current density is J(t) (J(t) = 
for t < 0) and the response j(t) is defined by 



j(t) = ^- j(t) = for t<0 (2) 

Then the complex conductivity cr(u) and the response j(t) = J(t)/E are 
related through 



oo 

-,•) j<'~'.,in<i/ (3) 





From (2), (3) one deduces that 
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Rea(u) = Rea(-w) 



(4) 



+00 

j(\t\) = - f e^Re^dw (5) 

7T J 



In (3) the integral over the time t runs over (t > 0) only due to causality 
(j(t) = for t < 0). This implies that the conductivity cr(u) is an analytical 
function of frequency u in the upper half of the complex plane. From the 
analyticity of a(cu) the Kramers-Kronig relations, that relate the real part 
and the imaginary part of the conductivity, can be deduced 



00 



lma(cu) = - PP / K -^du (6) 

7T J UJ — U 



-00 



ReaH = -ipP 7^Md. (7) 

7T J UJ — U 



-co 



where PP means the principal part of the integral. This implies also the fol- 
lowing spectral decomposition for z in the upper half complex plane: 



% r tie a(uj ) , , 

a(z) = - / Vd U ' 8 

7T J Z — UJ 1 

—00 

Finally we recall that the conductivity obeys sum rules. For example the re- 
sponse j(t = 0) is independent of the quantum character of electrons. It de- 
pends only on their concentration n, mass m and charge e through 



J(t = 0) = — (9) 
m 

Combining with (4), (5) one get 



+00 2 

Re a(u) du = (10) 





/ 



2.2 Relation between low frequency conductivity and quantum diffusion 



The quantum diffusion of states having an energy E is defined as AX 2 (E,t): 
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AX\E,t) = ([X(t)-X(0)f) E 



(11) 



where < A >e means an average of the diagonal elements of the operator A 
over all states with energy E. X(t) is the position operator along the axis x 
expressed in the Heisenberg representation. The velocity operator is defined 
as V x (t) = dX(t)/ dt, its correlation function C(E,t) is defined as 

C(E,t) = (v x (t)V x (0) + V x (0)V x (t)) E = 2Re(v x (t)V x (0)) E (12) 
and is related to quantum diffusion [6] through 



^AX 2 (E,t) = Jc(E,t')dt'. (13) 



As shown in [6], the real part of the low frequency conductivity is related to 
quantum diffusion. Indeed from the Kubo-Greenwood formula the real part of 
the conductivity is given by 




fi—hu> 



where \x is the chemical potential. In (14) the Fermi-Dirac distribution func- 
tion is taken equal to its zero temperature value. This is valid provided that 
the electronic properties vary smoothly on the thermal energy scale kT. For 
finite temperature, the effect of the Fermi-Dirac distribution function on the 
transport properties has been studied in the literature [7,8,9,10]. But, the- 
ses analyzes could not explain the unconventional conduction of quasicrystals 
and related alloys (very high resistivity at low temperature, and conductivity 
that increases strongly when defects or temperature increases). Therefore in 
the following, the Fermi-Dirac distribution function is taken equal to its zero 
temperature value. But the effect of defects and temperature (scattering by 
phonons...) on the diffusivity is taken into account via the relaxation time 
approximation (section 2.3). The central quantity F(E,u>) is given by 



AIT fiC I \ 

F(E, u) = —^-Tr(5(E - H)VJ(E + fuv - H)V X ) (15) 

where Q is the volume of the system and Tr means the Trace of an operator. 
Expressing the operator S(E — H) as the Fourier transform of the evolution 
operator e~ lHt one shows that 
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2nE - aJ] / <!/ v"' (YAt)VAO)) E (10) 



e 2 n(E) 
and 



2F(E-ftw,uj) 
e 2 n(E) 



/ dt e™(v x (0)V x (t)) E (17) 



where n(E) is density of states per unit volume (summed over up and down 
spins which are assumed to have the same transport properties here). Then 
one finds 



2Rea(E,w) = F(E,u}) + F(E-hu},uj) (18) 
where 



a{E,uj)=e 2 ^- Je iu,t C(E,t)dt (19) 
o 



Let us note that the function a(E, uj) is analytical in the upper half of the 
complex plane. For large uj : cr(E, u) oc 1/uj and the Kramers- Kronig relations 
are valid. Finally the usual sum rule is valid 



j Re9{E:Ul)dw= -J^El C{E:t = 0) = ^l (20) 



7re 2 n(E) 7ie 2 n 
nr a [Hi, uj) (±uj = " " 

o 

where m* is the effective mass and n the density of conduction electrons. 



If the variation of F(E, uj) with energy is small in the interval [E F —fkj, E F +hw] 
of values of E, one deduces from the previous set of equations that 



/ 7-1 \ 

Re a(u) ~ e 2 !^p> Re J e iu,t C{E F , t) dt (21) 

o 

(21) is valid at sufficiently small values of uj. In particular at zero frequency 
the dc conductivity is given by the Einstein relation 



o-(O) = e 2 n(E F )D(E F ) (22) 

with 
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D(E F ) = Um ^AX 2 (E F , t) (23) 



Finally there is a simple relation between the velocity correlation function at 
the Fermi energy and the impulse response j(t). Indeed comparing (21) and 
(3) one deduces the following equivalence at large time 



j(t) c e 2 ^C(E F ,t) (24) 
2.3 Relaxation time approximation (RTA) 



Within the relaxation time approximation one assumes that the response cur- 
rents respectively with disorder j(t) and without disorder j (t) are related 
through 

m=Jo(t)e-^ (25) 

where r is the relaxation time. So the relaxation time approximation (RTA) 
allows to treat the effect of disorder on quantum diffusion and conductivity. 
We give the relations satisfied by conductivity and quantum diffusion in this 
approximation. The conditions of validity of the RTA are discussed in part 3. 

Using (3), and within the RTA, the conductivity with disorder a(u,r) and 
without disorder ao(z) are related by 

a(u, t) = a (u + (26) 

The real part of conductivities with defects Re a(u, r) and without defects 
Re o"o(cu) are related simply. Using (3), it is straightforward to get 

Re a (u, r = — / £1-^- du' 27 

TXT J [UJ - UJ'Y + A 



which allows to compute the real part of the conductivity with defects. 

We discuss now the RTA from the point of view of quantum diffusion. In 
all cases we consider that the influence of disorder is much stronger on the 
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quantum diffusion than on the density of states. We thus neglect the variation 
of n(E) with disorder. From (24) one deduces that, for not too large disorder 
i.e. for sufficiently large relaxation time r the RTA is equivalent to 



C(E, t) = C (E, t) e-l*l/ T (28) 

where C(E,t) and Co(E,t) are respectively the velocity correlation functions 
with and without disorder. After equation (13) one deduces that the long time 
propagation is diffusive with a diffusion coefficient defined as 



+oo 

D(E) = -J C (E,t)e-W T dt (29) 
o 

which is equivalent to 

D(E) = ~AX 2 (E,t) if t»r (30) 



At zero frequency the diffusivity can be written in the useful form 



D(E*,r) = ^H (31) 



Using the t = conditions AX 2 (E,t = 0) = and f t AX 2 (E,t = 0) = and 
performing two integrations by part one get 



+oo 

J AX 2 {E F ,t) e- t/T dt 

L\E F ,r) = ° — = (AX\E F ,t)) T (32) 

J e^' T dt 



where {..J is a time average on a time scale r. AX (E,t) is the spreading of 
states of energy E, in the perfect system i.e. without disorder. 

More generally at low frequency, using (21) one can define a frequency depen- 
dent diffusivity D(E f ,uj) such that 

Re a(u) ~ e 2 n{E F ) D{E F ,u) (33) 
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and: 



1 °° 

D(Ef,w) = - Re J e iuJt C(E F ,t)dt (34) 


within the RTA (34) writes 



+00 



D(E F , u) = ^ Re J e (i "- 1/T)t C (£, f) df (35) 



It can be convenient to use the equivalent form which expresses the frequency 
dependent diffusivity D(E F ,cu) in terms of the quantum diffusion without 
disorder AX$(E,t): 



D(E F , u) = - Re H - - icuj y e^- 1 ^* AX 2 (£, t) dt I (36) 



^.^ Application to periodic Hamiltonians 



In this section we analyze the quantum diffusion in a perfect crystal, then we 
derive formulas for the low frequency conductivity within the RTA. 

Quantum diffusion 

Due to the Bloch theorem an eigenstate of a periodic Hamiltonian is defined 
by its wave vector k and by its band index n. The diagonal element of the 
velocity correlation operator on a state nk can be decomposed as follows: 

C(nk,t) = 2£ \v n>m {k)\ 2 cos (jE n (h) - E m (k)) . (37) 



where V is the velocity operator (in the chosen direction X) and V n>m (k) is 

— * — * 

the matrix element between states nk and mk. 

At small time, t <ti H/W where W is a typical bandwidth, using (13), the 
quantum diffusion is always ballistic with: 



AX 2 = Vlfi if t « A (38) 
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and 



V t 2 0t = (^2\(nk\V x \mk) 



(39) 



E„—E-p 



But in general the relevant time scale for electronic conductivity, which is 
the scattering time, is much larger than h/W. The following decomposition is 
important. 

Using (13) one shows that quite generally the quantum diffusion AX 2 (E, t) 
can be decomposed in a ballistic contribution and a bounded part: 



AX 2 (E,t) = V B (E)H 2 + AX 2 B (E,t). 



(40) 



The ballistic term VB{E) 2 t 2 is due to the diagonal elements of the velocity 
correlation function (intraband contribution) , whereas AX^ B (E,t) is due to 
the off diagonal terms of the velocity correlation function (interband contri- 
bution). One has 



Vb(E) = ( (nk\V x \nk) 



E n —E 



(41) 



and 



AX 2 B (E,t)<2AX 2 ooNB (E) 



(42) 



An important relation exists between AX^ B (E) and the square length of 
the unit cell in the chosen direction L 2 namely: 



(43) 



Indeed one shows easily that 



l-cos((E n -E m ){ 



AX 2 B (E F ,t) = 2h 2 ( £ , 

\m(m^n) V^™ Cjm > 



(nk\V x \mk) 



(44) 



E„—Ep 



the above expression is bounded by considering that all cosines are equal to 
— 1 and we define AX ( ^ )NB (£') as 
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\{nk\V x \mk) 



(E n — E m y 



(45) 



E n —E 



From the definition of the velocity operator 



X,H 



(46) 



and by considering the matrix elements of the position operator: 



(nk\X\mk) — J w* ^(f)xu m ^(r) dr 

Cell 



(47) 



one get 



AX^ NB (£) = 2( £ \(nk\X\mk) 



(48) 



E n —E 



Let us define an operator X~ that is constant in each unit cell and gives 
the average position of each unit cell along the chosen X direction. Since this 
operator is constant in each unit cell the orthogonality of mk and nk with 
{n 7^ m), implies: 



J u* n ^(r)x u m ^(r) dr = for n ^ m 



(49) 



Cell 



Thus the operator (X — X ) has the same off diagonal elements as X between 
mk and nk with {n ^ m) 



(nk\X\mk) = (nk\X - X-\mk) = J u* nR (f)(x - x~)u ml ;{r) dr (50) 



Cell 



— * 

The operator X P = (X — X~) has also well defined diagonal elements (nk\X — 
X~\nk) in the basis of Bloch states contrary to the operator X which has not 
well defined diagonal elements in this basis. 

Thus one can write 
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AXl NB (E) = 2 ( £ \(nk\X P \mk)\ \ (51) 



and 



AX^ NB (E) = ( (nk\X 2 P \nk) - (nk\X P \nk) 2 j (52) 

\ / E n =Ep 



AX^ NB (E) = ( (nk\(X P - (nk\X P \nk)) 2 \nk) j (53) 

The above expression depends only on the density probability of the Bloch 
wavefunction n ^(r)\ 2 . A bound of the above expression is easily established 

and according to the equation announced in (43) writes AX^ NB (£) < {^j . 
Low frequency conductivity in the RTA 

Let us consider now the low frequency conductivity and diffusivity of a crystal 
within the RTA. Using (40) and (36) one has 



D(E,r) = V B (Eyr+ LU ^ T) (54) 

It 

with 

V B (E) 2 = (\(nk\V x \nk) 2 ) (55) 



E n — E 



and 



J AX 2 B (t)e- t/T dt 

L 2 NB (E, r) = (56) 

/ e-*/ r dt 
o 

At low frequency, neglecting any variation of density of states with energy on 
a scale ftw, one can still write 

Re a(E F , u) = e 2 n{E F )D{E F , u) (57) 
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where the frequency dependent diffusivity can here also be decomposed in a 
Boltzmann and a Non Boltzmann contributions: 



D(E F ,u) = D B (E F ,u) + D NB (E F ,u) (58) 

with the Boltzmann contribution Db(E f ,u) and the Non Boltzmann contri- 
bution D^ B (E F ,u) given by 



D b (E f ,lu) = ^Re |(i -iuj J ^- l/T)t ^Xl{E F ,t)dt^ (59) 



2 + °° 



or 



D B (E ei u) = T ^- (60) 



and 



1 I 1 2 +00 I 

Z} NB (£ F ,^) = -Re <M--iu;) / e^-^AX^B^,*)^! (61) 
or equivalently 

D NB (E F ,u;) = ^Re | J e^' 1 ^ C NB (E F , t) d*J (62) 

where Cnb(£f,0 is the Non Boltzmann contribution to the velocity correla- 
tion function i.e. the off diagonal contributions n^min (37). 



#.5 Application to quasiperiodic Hamiltonians 



In this section we describe briefly the nature of states in a perfect quasiperi- 
odic system and the associated anomalous diffusion mode. Then we derive 
expressions of the conductivity within the RTA. 

Critical eigenstates in quasiperiodic Hamiltonians 

A first point is to define what a quasiperiodic Hamiltonian is [11,14]. Here we 
refer to the construction by the method of the acceptance zone in a crystal 
of higher dimension with a cut by a space with irrational slopes. We assume 
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also that the basic pattern, in the crystal of higher dimension, does not con- 
nect different unit cells as it is the case for some incommensurate phases. An 
important consequence is that a local environment of size L must be repeated 
exactly within a distance of order L. 

Consider a general one dimensional tight-binding Hamiltonian H of the form 



# = £^>01 + £^)<;| (63) 

where the first part corresponds to hopping between nearest neighbors 
and the second part to on-site energies e^. For the Fibonacci chain: 



ei = V[iip\ (64) 
V(x) =V for - up < x < uj 3 (65) 
V(x) =V 1 for - y? 3 < x < y? 2 (66) 

where (p is the Golden mean. For this Hamiltonian all states are critical. A 
critical state is intermediate between spatially extended and exponentially lo- 
calized. Its envelopp presents an algebraic decay with distance [12,13,14,15,16]. 

It is interesting to compare the eigenstates of the Fibonacci chain with those 
of two classical models: the ID Anderson Hamiltonian and the Harper Hamil- 
tonian. For the Anderson Hamiltonian the on-site energies e» are chosen ran- 
domly between two values [— W, +W]. In that case all states are exponentially 
localized with distance from some point. The localization length C depends 
on the energy and on disorder (through the parameter W/t) and decreases 
when W/t increases. For the Harper Hamiltonian the on-site energies are of 
the form 



€i = A COs(2"7TCc«) (67) 

For uj rational the structure is periodic and the eigenstates are Bloch states. 
For uj irrational the structure is aperiodic. Depending on the parameter A the 
states are either extended (A < 2), critical (A = 2) or exponentially localized 
(A>2). 

A remarkable difference between the Harper and the Fibonacci model is that 
in the second model the states are critical for any values of the parameters. It 
is clear that the exact repetition of a given environment has a strong influence 
on the long range correlations of eigenstates. Indeed in the Fibonacci model 
a given sequence of length L is repeated exactly within a distance of order L 
whereas in the Harper model (or in the Anderson model) this exact repetition 
does not exist. 
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Fig. 1. Critical eigenstate in a decagonal model. The wavefunction density prob- 
ability is small on white spots and strong on black spots. From T. Fujiwara [17]. 



A major question is whether the critical states persist in higher dimension. 
Although there could be exceptions it seems that in general states are critical 
in 2 or 3 dimensions. There have been numerical studies of spectra in more 
than one dimension and figure 1 presents a state calculated for a large 2- 
dimensional model by T. Fujiwara et al. [17]. 

An argument proposed by C. Sire [18] is the following. If a wavefunction \&^ 
is mainly localized in a region of extension Lq then it should live also on 
any similar environment. After Conway theorem, identical environments are 
located at a distance 2L at most. Introducing a tunneling factor z between 
the two local environments we write 



This qualitative argument points out the importance of the exact repetition of 
an environment within a distance comparable to the size of this environment, 
which is typical of a quasicrystal. 

Band scaling and anomalous diffusion 



*2L = zm Lo and # 2 n Lo = z n y Lo 



(68) 



Introducing L 



2 n L , we can write equivalently 




(69) 
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The critical states are associated to a fractal energy spectrum in quasicrys- 
tal. The nature of this spectrum can be understood by considering a series 
of approximants with unit cell of increasing size L. When going from one ap- 
proximant to the following, the bands are broken into smaller pieces which 
gives a fractal spectrum for the quasiperiodic system. 

Typically the width AE of a band varies like 

AE ~ — (70) 

where L is the length of the unit cell of the system. The exponent 7 is greater 
than one due to the effect of the quasiperiodic potential. 7 depends on the 
energy of the wavepacket and of course on the parameters of the Hamiltonian. 
This scaling law reflects also the critical character of eigenstates and their 
algebraic decrease with distance. 

The band scaling has a direct consequence for the propagation of waves in 
quasiperiodic media. The spreading L(t) of a wavepacket is neither ballistic 
(i.e. proportional to time t) as in crystals nor diffusive (i.e. proportional to 
the square root of time) as in disordered metals. In general it follows a power 
law and one can write after a time t: 



L(t) ~ At 13 (71) 

/3 depends on the hamiltonian and on the energy of the wavepacket [19]. Indeed 
for an approximant with unit cell of size L the characteristic velocity v(L) is 
given by 



, T . 1 dE(k) LAE B 
v(L) oc y- oc — — oc — — 72 



A; is a wave vector in reciprocal space. When the magnitude of Ah is of order 
1/L (i.e. the size of the Brillouin zone) the magnitude of the energy variation 
is of order AE (i.e. the width of a band). The last equality in (72) is then 
obtained through equation (71). Since 7 > 1 the typical velocity tends to zero 
when L increases. 

Using the relation between the spreading of the wavepacket and the velocity 
at length scale L: 



— — = v[L) oc - — - (73) 
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One integrates straightforwardly the above differential equation and obtains: 



Lit) ~ At 13 with = -<l 

7 



(74) 



Thus in a quasiperiodic system the asymptotic diffusion law is at sufficiently 
large t: 



but let us recall that at small time the propagation is necessarily ballistic. 
Indeed after (13) and (12) one has 



There is numerical evidence of anomalous diffusion in 1- or 2-dimensional 
systems (see figures 2, 3). It should be noted that these law are numerically 
approximate because in practice there is always fluctuations. As a rule, one 
expects that the fluctuations are less important when the dimensionality of 
the system increases. 

Associated to this anomalous diffusion one observes also an anomalous trans- 
mission of waves through a finite part of a quasiperiodic tiling of length L. 
The transmission coefficient varies like a power law of the length L. There is 
also a power law variation of the resistance of a stripe of length L. 

It is interesting to consider the propagation law AXq (L, t) in an approximant 
of unit cell size L, in the limit of large L, and to relate it to the propagation 
in the quasicrystal AXg(t). 

In the periodic approximant of unit cell size L after equation (40): 



AA 2 (t) ~ At 20 



(75) 




(76) 



AX 2 (L,t) 



V(L) 2 t 2 + AX 2 B (L,t) 



(77) 



or equivalently 



AX 2 



NB 



(L,t) 



AX 2 



(L,t)-V(L) 2 t 2 



(78) 



Since the velocity V(L) tends to zero at large L one gets 
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Fig. 2. Behaviour of Rq(E, t) = AX 2 (E,t)+AY 2 (E,t)+AZ 2 (E,t) for three energies 
in the 2- (left) and the 3-dimensional (right) generalized Rauzy tiling. From [20]. 
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Fig. 3. Normalized density of states (lower curves) and diffusion exponents (upper 
curve) of the 2- (left) and the 3-dimensional (right) generalized Rauzy tiling. From 
[20]. 



lim AX 2 B (L,t) = lim AX 2 (L,t) = AX 2 {t) 



(79) 



Thus at a given time t the Non Boltzmann contribution to the spreading 
AX^ B (L,t), for an approximant of unit cell size L, tends to the anomalous 
diffusion law of quasicrystals AX^t) in the limit of large L. 



Low frequency conductivity in the RTA 
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According to the anomalous law derived above and exemplified by numerical 
simulations we assume that at sufficiently large time t — > oo the diffusion law 
for a perfect quasicrystals is 

AX (E,t) 2 ~ A(E)t 2f3{E) (80) 

Let us recall that at small time t the diffusion law is ballistic (76). According to 
the general expression relating frequency dependent diffusivity and quantum 
diffusion (36) one has in the RTA: 



2 +°° 



Rea(E F ,u) = ^f^-Re | " ^) / e^' 1 ^ AX 2 (E, t) dtj (81) 

If the asymptotic law (80) is applicable on the time scales r and 1/uj (which 
means a limit of small disorder and small frequency) then one can write 

Re a(E F , u) = Re a(E F , u) (82) 

with 

^,.) = ^r(2, S + l)(^f' (83) 

2 VI — IUJT J 



where V is the Euler Gamma function. The formula (82), (83) is a generalized 
Drude formula. Indeed if one considers the case (3 — 1 then one recovers 
exactly the Drude formula. 

The generalized Drude formula is established for a given (3 in the limit of 
infinite scattering time r — > oo and low frequency uo — > 0. Thus this formula 
cannot be applied for a given r in the limit /3 — > 0. In order to treat the case 
of small (3, for a fixed r, we proceed in the following way. We start from the 
low frequency form (84) for the perfect system and treat the effect of disorder 
within the RTA. For the perfect system we assume the conductivity satisfies: 

/, ,\ 1-2/3 

Re<7 (c<j) = <7 ( — j for < u>i (84) 



After (3), (25) the conductivity a(uo,r) of the system with defects satisfies: 
a(u, t) ~ a + (85) 
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thus cr(u, t) is determined once gq{z) is known for z = cu + i/r. In the limit of 
low frequency and long scattering time we need a (z) for small values of \z\. 

Consider the two functions that are analytical in the upper half of the complex 
plane: 

a / 7. \ !- 2 ^ 



<7 x (z) = &(z) for x > (87) 



In (86), we choose the following determination for the power of a complex 
number. For a complex number X = \X\e' ld with — n < 9 < n we define 
X a = \X\ a e iad . From (84,86), Rea (uj) = Rea(cu) for \u\ < ooi. 

For a function f(z) which is analytical in the upper half of the complex plane, 
and which modulus tends to zero at large z one has 

7T J Z — UJ 

— oo 



The spectral representation (88) applies to both cr (z) because it is the con- 
ductivity of a real system and to a x (z) because of the factor ixj (ix + z) which 
guarantees that the function a x (z) is analytical in the upper half complex 
plane (x > 0) and decreases sufficiently quickly at large z. So 

ao(z) = o x {z) + % - + r ^^oM-^^M duj m 

TV J Z — UJ 

— oo 



Since Re<7o(o>) and Kecr x (uj) are pair functions of uo one can add together 
contributions of u and — u in (89). Using Recro(^) = Rea(w) for |u;| < uj\ and 
taking the limit x — > +oo one gets 

. . 2zi + r Rear)(cu) — Rea(cu) , . , 

„„(,) = „(,)-_/ J (i >_ (90) 



Developping the integral in a power series of z and keeping only the term 
proportionnal to z one gets 

/ \ ~t \ izoQ 2izaa ( ( z \ 3 \ 

^»='< :)+ ^— +o (U) (91) 



23 



a 



uj\ /' Re<xo(u;) 



O"0 



/ 



da; 



or 



(92) 



For a given /3 7^ the first term in the right member of (91) dominates in the 
limit z — 0. In this limit, one recovers, as expected, the Generalized Drude 
formula (82), (83) and [6]. 



Yet it is interesting to look at the limit of (3 — > for a given r or given 
frequency u> i.e. for a given value of z = uj + i/r. In this limit, not considered 
in [6], one gets from (91), (82) and (83): 



2u\z 



a - log ( — J 



(93) 



In (93) log(X) = log I A" I + i9 for X = \X\ e ie with -tt < 9 < tt. Introducing 
A = 2o"o/tt^i the dissipative part of conductivity, up to terms of order (z/cui) 3 , 
is given by 



Re a ~ 



a + log 



0J\T 



+ wrArctg(c<jr) 



(94) 



Finally we note that the quantum diffusion associated to the conductivity 
law (94) is logarithmic. Indeed the zero frequency diffusivity is given by D oc 
j (a + log(cjir)) and after (31) the mean free path is L 2 (t) oc (a + log(c<Jir)). 
Physically the mean free path represents the diffusion length between two 
scattering events in the perfect system. 



3 Anomalous quantum diffusion and conductivity in periodic and 
quasiperiodic systems 



In this part 3 we discuss first the limitations of the RTA and the possibility 
of a metal-insulator transition. Then we focus on the physical interpretation 
and consequences of the formulas derived in part 2. We apply this to a general 
theory of low frequency conductivity and metal-insulator transition in crystals 
and quasicrystals. 

Let us recall first the main formulas derived in part 2 to treat low frequency 
conductivity. 

In the impulse response formalism one considers the response of a system to 
an impulse of electric field E(t) = E5(t) where 8(t) is the Dirac function. The 
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resulting current density is J(t) (J(t) = for t < 0) and the response j(t) is 
defined by j(t) = J(t)/E. The complex conductivity cr(u) and the response 
j(t) are related through 

oo 

<t(w) = J e iu}t j(t)dt (95) 



The quantum diffusion of states of energy E is measured by the quantity 
AX 2 (E,t): 



AX\E,t) = {[X(t)-X(0)]') E (96) 



where < A >e means an average of the diagonal elements of the operator A 
over all states with energy E. X(t) is the position operator along the x axis in 
the Heisenberg representation. The velocity operator V x (t) = dX(t)/ dt has a 
correlation function C(E,t) defined as 



C(E,t) = {V x (t)V x (0) + V x (0)V x (t)) E = 2Re{V x (t)V x (0)) E (97) 



and is related to quantum diffusion through 



d 

~dt 



AX 2 (E,t) = j C(E,t')dt'. (98) 



The low frequency conductivity satisfies 

Re a(u) ~ e 2 n{E F ) D(E F ,u) (99) 
with n(E) the total density of states and 

D(E F ,u) = -Re J e iw)t C{E F ,t)dt (100) 
o 

In the limit of large time j(t) and C(E F ,t) are related by 
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j(t)^e*^C(E F ,t) (101) 

where E F is the Fermi energy and n(E F ) the density of states at the Fermi 
energy (summed over spin up and down). 

Within the relaxation time approximation (RTA) one assumes that j(t) and 
C(E, t) with disorder are related to j (t) and Cq{E, t) without disorder through 



j(0=Jo(*)e-" l/T (102) 



C(E,t) =C {E,t)e-W T (103) 



Here the Fermi-Dirac distribution function is taken equal to its zero tempera- 
ture value. This is valid provided that the electronic properties vary smoothly 
on the thermal energy scale kT. But in the RTA, the effect of defects and 
temperature (scattering by phonons...) is taken into account through the re- 
laxation time r. The diffusivity is given by 



+oo 

D(E f ,uj) = -Re J e {i "- 1/T)t C (E F ,t)dt (104) 



It can be convenient to use the following equivalent form which expresses the 
frequency dependent diffusivity D(E F ,cu) in term of the quantum diffusion in 
the systemwithout disorder AX$(E, t) 



+oo 

D(E f ,uj) = -Re \(- -iu) [ e^' 1 ^ AX^(E F ,t)dt\ (105) 




o 



Note that (100), (105) are valid at sufficiently small frequency. The condition 
is that the conductivity varies only slowly with E F on the energy scale Huj. If 
this is not the case it is still possible to define an average of the conductivity 
over a range of values of the energy E F . In that case the above formulas are 
applicable to the conductivity averaged on a scale of Fermi energies greater 
than hu. 

At zero frequency the above formula (105) is exact within the RTA and can 
be written as 
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Fig. 4. Interference between a path (solid line) and its time reversal image (dashed 
line) for an electron which diffuses through a static disordered potential. The in- 
terference that occurs between two different paths after several scattering events 
cannot be described by the RTA. 

D= LH f; T) (106) 

with the mean free path L(E F , r) given by 



+00 

L 2 (E F , t) = - f Q~ t/T AX*(E F , t) dt (107) 



3.1 Validity of the RTA and Anderson transition 



Let us discuss now the conditions of validity of the RTA. We consider first 
the role of elastic scattering i.e. diffusion by a static potential. From the def- 
inition of the RTA it is clear that j(t) = j (t) exp (— \t\/r) is essentially zero 
beyond the the relation time t > r. However this is not the case in many disor- 
dered systems which present elastic scattering. For example in the case of free 
electrons scattered by static defects at random in a three dimensional system 
there are interferences between several scattering paths. This is represented 
in figure 4. As a consequence of these interferences the long time behavior of 
the response current j(t), in the absence of a magnetic field, is 



j(t) = -Ar 3/2 with A > (108) 



Thus interferences effects in the diffusive regime cannot be in general described 
properly by the RTA. Indeed in the RTA all quantum correlations are lost 
beyond the scattering time scale r. Figure 4 shows a type of interference that 
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occurs on a time scale greater than the scattering time r and thus cannot be 
described in the RTA. 

Ultimately these interferences can lead to a localization of the states, provided 
that the disorder is sufficiently strong. This Anderson localization is schemat- 
ically represented in figure 5. Of course a correct theory of the Anderson 
localization is also beyond the RTA. 

For a 3-dimensional system the importance of quantum interferences depends 
on the ratio between the dc-conductivity of the system cr dc and the Mott 
value o"Mott — 600(ficm) _1 /A where A is the mean-free path expressed in 
Angstroms. If R = o del® Mott 1 the effect of the quantum interferences on 
a dc is small. In that case the quantum interferences in the diffusive regime 
change only slightly the value of the conductivity and the RTA can describe 
the role of elastic scattering. If the ratio R is closer to one then the RTA 
cannot be used to describe the role of elastic scattering. 

The case of inelastic scattering is different. It is generally assumed that in- 
elastic scattering with a scattering time Tj n destroys quantum interferences 
on a time scale t > Ti n . So the relaxation time approximation is expected to 
be valid in the case of inelastic scattering. Yet a condition is that the system 
contains a sufficiently large number of states so that the electron can find a 
state to scatter into close to the scattering event. 

For this reason the hopping regime between localized states (either short range 
hopping or variable long range hopping) is not described by the relaxation time 
approximation. Indeed each hopping process requires an exchange of energy 
with phonons that provides to the electron the difference in energy between 
the initial and final localized states. 

General case according to the scaling theory 

We discuss now the case where there is both elastic scattering (characterized by 
r e i) and inelastic scattering scattering (characterized by Tj n ) and assume that 
r in > T e i. This often happens since Tj n diverges at low temperature whereas r e i 
is temperature independent at sufficiently low temperature. In the case where 
T %n < T e i the elastic scattering is expected to have a minor effect on transport, 
and one is back to the case of pure inelastic scattering. 

Let us define t(L) the time needed for a wavepacket to spread on a length 
scale L and D(L) = L 2 /2t(L) the diffusivity D(L) at length scale L. Let us 
define also L in which is such that r in = t(L in ). 

According to the scaling theory the diffusivity depends on the length scale 
due to the quantum interferences in the diffusive regime. One can distinguish 
three steps. 
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Fig. 5. Schematic band structure of an Anderson insulator, with increasing disorder 
from bottom to top. At zero temperature the system is metallic if the Fermi energy 
is in the extended states region. The system is insulating if the Fermi energy is in 
the localized states region. 

1) Define the diffusivity at the length scale of the elastic mean free path. This 
can be done by using the RTA with a relaxation time equal to the elastic 
scattering time. This is valid because the quantum interferences do not act on 
length scales smaller than the elastic mean free path. 

2) Consider the conductance g(L) of a cube with size L. 



g(L) = e 2 ^^y = e 2 n(E)D(L)L (109) 



and apply the scaling relation on length scales L such that L in > L > L e . 
dlog(fif) 



dlog(L) 



P(g) (no) 



beta(g) is represented schematically in figure 6. The relation (110) defines 
entirely the quantum diffusion and thus the velocity correlation function as a 
function of the length scale L or also as a function of time scale t(L). Thus it 
allows to compute C sca ii ng (E,t) where C sca u ng (E,t) is the velocity correlation 
function corresponding to the purely elastic scattering. 

3) One stops the renormalization procedure at the length scale L in = L(T in ), 
which means that the macroscopic diffusivity is simply D = D(L in ). Equiva- 
lently the role of inelastic scattering is simply to destroy the velocity correla- 
tion on the time scale Tj„ i.e. one has 
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Fig. 6. Schematic representation of the (3(g) function for systems of dimension d. 

C(E,t)=C scalmg (E,t)e- t ^" (111) 

The renormalization procedure determines entirely C sca i ing (E, t) and thus through 
(111) determines entirely the quantum diffusion of the system with inelastic 
scattering. 



3.2 Phenomenon of backs cattering 

Let us consider the response j(t) to an electric field E5(t) applied to a system. 
During the impulsion the dynamics is dominated by the effect of the electric 
field. This implies that j(0) is independent of the atomic forces applied to the 
electrons and is given by the classical response j(0) = ne 2 /m. In particular 
j(0) is positive. In the classical Drude model j(t) = j(0)exp (— t/r) tends to 
zero at larger times but is always positive. This is illustrated by figure 7. 

Quantum effects can lead to a counter intuitive situation where j(t) becomes 
negative in some time interval. This is the phenomenon of backscattering. This 
phenomenon occurs in a recurrent manner in the physics of quasicrystals and 
related complex intermetallics. So we give here an account of its characteristics 
and its consequences on conductivity. 

We illustrate first the phenomenon of backscattering in two cases concerning 
disordered systems: the weak-localization and the strong localization regimes. 

Inelastic scattering in the weak localization regime 
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Fig. 7. Response current j(t) within the Drude model. Without disorder (solid line) 
the response is constant j = ne 2 /m. With disorder (dashed line) the response decays 
exponentially on the time scale r. 
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Fig. 8. Response function j(t) in a disordered system with weak-localization effect, 
and thus backscattering at large time. Without inelastic scattering (solid line) and 
with inelastic scattering (dashed line). 

We assume that the elastic scattering time r e is much shorter than the inelastic 
scattering time r in . On the time scale r e < t < r in the interferences are 
insensitive to inelastic scattering and j(t) = —At^l 2 (for a 3-dimensional 
system), but on the time scale t > Ti n the interferences are destroyed. Thus a 
general expression of j(t) for t > T e \ is (figure 8): 



j{t) = -Ar z/2 e~ t/Tin with A > 
In that case using equation (95) one finds for the dc conductivity: 

2A^ 



(112) 



= Jj(t)dt = 



O"0 + 



(113) 



where a is the dc-conductivity in the absence of inelastic scattering. That is 
inelastic scattering tends to increase the conductivity. This paradoxal result 
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Fig. 9. Response function j(t) for an insulator. Beyond a characteristic time 
(T c ~ 2 here) the response jo(t) is essentially zero. Without inelastic scatter- 

oo 

ing (solid line) the conductivity is zero and the integral a dc = f jo(t)dt = 0. 

o 

With inelastic scattering (dashed lines) the system is conducting and the integral 

oo 

<Jdc = /jo(t)e-*/ r ™ dt >0. 



is directly related to the backscattering phenomena that occurs with static 
disorder. Indeed in the integral (113) with j(t) = —At~ ? 'l 2 e~ t l Tin , r in can be 
considered as a cut-off which supresses the negative contributions of j(t) = 

-At~ 3/2 &tt> T^. 

The role of a small frequency u> is comparable to that of inelastic scattering 
at a time scale r in ~ 1/lu. It acts also as a cut-off. For example in the case of 
the weak localization one has for r,„ — > oo: 



Re a(u) = o- + Ay/2iru (114) 

This means that Re a(u, r) increases with frequency at small frequencies. This 
is at the opposite of the standard behavior of metals, where the low frequency 
conductivity is dominated by the Drude peak. 

The role of frequency can also be deduced from the relation (27) between 
the real parts of conductivity with scattering and without scattering. Indeed 
this relation (27) shows that is the average of Reao(c^) by a Lorentzian 
of width 1/r centered at zero frequency. a dc can increase with 1/r only if 
Rea (u) increases with uj. 

Inelastic scattering for localized states : the Thouless regime 

Let us consider a system that it is an insulator. In that case the zero frequency 
conductivity is zero: 
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0)= Jj(t)dt = (115) 



Since the integral is zero, and since the small time response j(t) > it means 
that there are times interval where j(t) < 0. Thus in an insulating system the 
backscattering phenomenon necessarily occurs. For Anderson localization one 
expects that j(t) > at small times and j(t) < at large time. Beyond a 
characteristic time T c , the response j(t) is essentially zero (see figure 9). 

Provided that the RTA is applicable one get 



oo 

a(uj = 0, r) = J j (t) e~ t/T dt = (116) 



If r > T c then the exponential differs significantly from 1 only in the large 
time region t > T c . In that case one can replace exp(— t/r) by 1 — t/r in the 
range (t < T c ) where j(t) ^ 0. This gives 



1 °° 

a(u = 0,T) = --Jj o (t)tdt = (117) 



Although J j (t)dt = the integral J j (t)tdt in (117) is negative because 

o o 
jo(t) is negative at large times and positive at small times. 

The regime described here correspond to the Thouless regime, that can be 
described by analyzing the quantum diffusion. In the Thouless regime, the 
physical picture is that of electrons spreading during a time r in between two 
inelastic scattering events and then loosing completely their phase at each 
inelastic scattering event. According to the Thouless scenario, in the limit of 
large inelastic scattering time Tj„, the spreading of the electron wavefunction 
between two inelastic scattering events is bounded by the localization length 
£e- Since the electron must wait till the next inelastic scattering event to loose 
phase memory and spread again the diffusivity is given by 



£2 

DThouless(E) = (118) 



' in 



Within the RTA the diffusivity which is the square of the spreading during r i? 
divided by r in is given by 
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D (E) = L2(E ' T ^ (119) 
with the mean free path L 2 (E,T in ) given by (107) which is equivalent to 



L ' {E * ) = £ iE-WH*,*.)* Xi " <120) 

with 



X EiE , = (E\X\E') (121) 



At large inelastic scattering time r in one gets L 2 (E) = lim^^oo L 2 (E, r h 
which can be written as 



L 2 {E) = (E\(X - (E\X\E)) 2 \E) = £| (122) 

where £g is the localization length of the state of energy E. After (119) the 
diffusivity is D(E) = £ E /2ri n , in agreement with the argument of Thouless. 

One thus recovers the typical dependence on the scattering time r i.e. a^c oc 
1/r that was deduced from the analysis of the response j(t) and the backscat- 
tering. 



3.3 Anomalous quantum diffusion and conductivity of periodic systems 



Conductivity within the RTA 

The semi-classical theory of conduction in crystals is based on the concept of 
a charge carrier wave-packet propagating at a velocity V = (1/ h)dE n (k) / dk, 
where u E n (k)" is the dispersion relation for band n and wavevector k. The 
validity of the wave-packet concept requires that the extension L wp of the 
wave-packet of the charge carrier is smaller than the distance Vt of traveling 
on the time scale r. On the contrary, if Vt < L wp , the semi-classical model 
breaks down. The quantum formalism presented here allows to treat on the 
same footing the standard regime where the semi-classical approach is valid 
and the small time regime Vt < L wp . As shown in part 2 the spreading of 
states with energy E in crystals is given by 



AX 2 (E F , t) = V 2 t 2 + AX 2 B (E F , t) (123) 
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The first term in the right hand side of (123), V 2 t 2 , corresponds to the Boltz- 
mann contribution. This term dominates at large times and describes the in- 
tercell ballistic propagation of wavepackets on long time scale in crystals. The 
physical origin of this term is the coupling between successive unit cell that al- 
lows the electron to travel in the whole crystal. The second term AX^ B (E F , t) 
is the Non-Boltzmann contribution. It describes the intracell spreading of the 
electron. Indeed as shown in part 2 this spreading is bounded by a term of the 
order of the square of the unit cell size. In a standard crystal the Boltzmann 
term dominates at the time scale relevant for transport, i.e. the scattering time 
due to disorder In approximants of quasicrystals the Non Boltzmann term can 
dominate. 

Let us anticipate on the ab-initio calculations which show that AX^ B (E F ,t) 
is nearly constant AX^ B (E F ,t) ~ AX^ B except at very small times t <C r. 
Then the equation (105) leads to 



D(E F , u>) * + (124) 

Thus the frequency dependent diffusivity D(E F , uj) is the sum of a Drude like 
contribution (first term on the right hand side of (124)) and a contribution 
independent of frequency which increases with disorder. As we show from ab- 
initio calculations (see below) the Drude like contribution can be small in some 
periodic approximants of quasicrystals. This explains why in these systems 
the optical conductivity presents no Drude peak and why the dc- conductivity 
increases with disorder. 

We define r* as the time for which the Boltzmann and Non Boltzmann con- 
tributions to the dc- diffusivity are equal. Thus: 



r* = ^ (125) 



In that case the dc-diffusivity can be represented as in figure 10. The ac- 
diffusivity is represented in figure 11. 

Backs cattering 

The Non-Boltzmann contribution to the diffusivity D NB = AX^ B /(2r) is 
formally similar to the Thouless expression for localized states D Thou i ess = 
£b/(2t). Indeed in the Thouless regime the quantum diffusion between two 
inelastic scattering events is limited by the localization length For crystals 
it is the Non-Boltzmann contribution to quantum diffusion which tends to 
AX NB . Let us recall that AX NB is itself limited by a term of the order of the 
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Fig. 10. Variation of the normalized zero frequency diffusivity D/V 2 t* with scat- 
tering time for the expression (124). For r = r* the Boltzmann and non-Boltzmann 
contributions are equal. 












* T= T 
" T=T/3 
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Fig. 11. Variation of the normalized frequency dependent diffusivity D/D^b- For 
t > t* the transport is dominated by the Boltzmann term and the Drude peak is 
well defined. For r < r* the transport is dominated by the Non-Boltzmann term 
and the Drude peak is absent. 



unit cell size. 



As discussed previously (see section 3.2) the increase of conductivity with 
disorder is a direct consequence of the backscattering. Indeed we shall find in 
part 4 that ab-initio calculations prove the existence of backscattering. 

Finally we emphasize an important difference with the Thouless regime. In 
the Thouless regime it is the inelastic scattering that destroys the localization 
produced by the elastic scattering. Here, provided that the RTA is valid, it 
is either the elastic or the inelastic scattering that destroys the localization 
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induced by the periodic potential. 
Metal-Insulator transition 



Let us discuss now the role of quantum interferences according to the scaling 
theory of localization. As explained in part 3.1, the central quantity is the 
conductance of a cube with a size equal to the elastic mean free path L(E F , r). 



g ~ e 2 n(E F )D(E F , t)L{E f , t) (126) 



Since the quantum interferences effect have not the possibility to operate at 
smaller length scale than L(E F , r) then this quantity can be computed with 
the RTA according to (107). 

We still assume that AX NB (E F ,t) is nearly constant equal to AX NB except 
at the smallest time (see below the ab-inito results on an approximant of 
a-AlMnSi). 

Thus the typical propagation length L(E F , r) on a time scale r, i.e. the mean- 
free path, is 



L 2 (E F , r) = AX 2 B + 2V 2 r 2 (127) 

Let us introduce g , which is characteristic of the perfect crystal and is defined 
by 



g = e 2 n(E F )AX 2 B V (128) 
Let us introduce an adimensional value f of the scattering time r defined by 



f = ^t = ^? (129) 



Let us define also the function f(x): 

f{x) = {h +x ) ^ 1+2x2) (130) 

Then one has 
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Fig. 12. Metal-Insulator phase diagram as a function of the two parameters <?o/<?c 
and 1/f = ^1^. The insert represents the limit of a normal metal i.e. for fixed r 
and V the limit of a small AXnb- After (128) and (129) this limit is in the region 
of the phase diagram at small go/g c and small 1/f. 



9 = 9of(f) 



(131) 



After the scaling theory a three dimensional system is insulating (metallic, 
respectively) if g < g c (resp. g > g c ) where g c is the value of the universal 
critical conductance in the scaling theory. Using g = gof(f) it is equivalent 
to say that the system is insulator if go/g c < l/f(f) and metallic if go/g c > 
l//(f). This is illustrated in figure 12. Note that go/g c is characteristic of 
the perfect crystal whereas 1/f measures the scattering rate 1/r in units of 
V/AX NB . 

A first remarkable property of this phase diagram is that if go > Rg c with 
R = (2/3) 3 / 2 ~ 0.54 then the system is always metallic whatever the value 
of the scattering rate (phase (a) in figure 13). This is not the case for normal 
metals that always become insulating at sufficiently small scattering time r 
(i.e. at sufficiently large disorder). 

If g < Rg c the system is metallic at large and small scattering rates and 
insulator in an intermediate zone (phase (b) in figure 13). This means that 
if the system is in the large 1/f metallic region it will become insulating by 
decreasing 1/f that is by decreasing disorder! This is just the opposite of the 
standard conditions for the occurrence of the Anderson localization transition. 
The other insulator-metal transition is normal in the sense that the metallic 
state is obtained by decreasing disorder. 

Note that the case of a normal metal corresponds to the limit AX NB — > 0. In 
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Fig. 13. Two types of systems can exist. Phases of type (a) are always metallic 
whatever the value of the scattering time r. Phases of type (b) can be either metallic 
or insulating depending on the value of the scattering time. For phases of type (b) the 
Metal-Insulator transition that occurs at the highest value of 1/f is unconventional 
since the metallic state is obtained by increasing 1/f i.e. by increasing disorder. 

that case one uses the asymptotic form of the function /(f) for large f namely 
/(f) ~ \f2f 2 . One then recovers the standard criterion for free- like electrons. 

3.4 Anomalous quantum diffusion and conductivity of quasiperiodic systems 

Conductivity within the RTA 

A generalized Drude formula (82,83) for the low frequency conductivity is 
derived in part 2.5. Yet it is interesting to derive it from simple physical argu- 
ments. One notes first that at zero frequency the dependence on the scattering 
time is easy to establish. Indeed the diffusivity is 

D(ft,r) = ^%2 (132) 
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Fig. 14. Dissipative part of the conductivity R(u) normalized by the zero frequency 
value (R(co) = Kea(uj)/a(uJ = 0)), as a function of lot (frequency normalized by 
1/t). The left panel shows that the conductivity increases with frequency when 
(3 < 0.5. The right panel shows that there is a low frequency peak if (3 > 0.5. The 
case (3 = 0.5 is an intermediate case for which the conductivity is independent of 
frequency. 

where L 2 (E F , r) is the mean free path i.e. a typical distance of propagation in 
the perfect structure during the scattering time r. Assuming that in the perfect 
quasiperiodic structure the spreading of a wavepacket is AX 2 (t) ~ At 213 one 
gets 



L 2 (Ep, r) ~ At 213 



(133) 



and 



D(t) 



L 2 (E f ,t) _A 



2t 



2 



2/3-1 



(134) 



One notes also that the conductivity depends on scattering time and frequency 
only through the combination r/(l — iuor). This stems directly from the RTA 
formula (105) which is expressed as a Fourier-Laplace integral with (^ — iuj) = 
(1 — iujt)/t. Thus the frequency dependent diffusivity is 



D(u,r) ~ ^ R e ( , T ) 
v ; 2 \1-iuot) 



\ 2/3-1 



(135) 



Except for a numerical factor T(2f3 + 1) this formula is equivalent to the 
Generalized Drude formula (82), (83) which we recall here: 
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Re a(Ep, to) = Re a(E F , uj) 



(136) 



with 

ift.l^iW+Hm''" (137) 
2 VI — iujt J 

The behavior of the conductivity depends on the value of (3 compared to 
0.5. The frequency dependence is represented in figure 14. If f3 > 0.5 the 
behavior is similar to that of a normal metal. The dc-conductivity decreases 
when disorder increases and the low frequency conductivity presents a peak at 
low frequency, somewhat similar to the Drude peak. If /3 < 0.5 the behavior 
is not that of a metal. In the absence of disorder the system is insulating, and 
the dc-conductivity increases when disorder increases [24]. The real part of 
the conductivity increases when frequency increases. Instead of a Drude peak 
there is a dip. 

One notes also that even in the absence of scattering i.e. for r — > oo the 
real part of the conductivity is non zero in the limit of small frequency. This 
means that there is absorption of electromagnetic energy by the system. This 
is not the case in a normal metallic crystal. Here the absorption of energy can 
be understood by considering approximants with large unit cell. In perfect 
approximants the absorption of energy is made through interband transitions. 
For a given frequency the absorption of energy by interband transition can 
occur with sufficiently large unit cell because the bands become very narrow 
and very close in energy allowing for interband transitions. 

Backs cattering 

For (3 < 0.5 the behavior of conductivity with frequency and disorder is not 
that of a metal. This can be attributed to the phenomenon of backscatter- 
ing. Indeed after the relation (98) between quantum diffusion and velocity 
correlation one gets 

C (E, t) = d2AX ^ E,t) * Aim ~ (138) 

if the quantum diffusion law is AXq (E,t) ~ At 213 in the perfect quasiperiodic 
system. From (138) it appears that the velocity correlation function is neg- 
ative at large time if /3 < 0.5. The response j(t) is then also negative. The 
phenomenon of backscattering then implies that the conductivity increases 
with disorder or with frequency, as discussed for weak-localization regime or 
for the Thouless regime (see figure 15). 
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Fig. 15. Typical behaviour of the response j(t) for a quasiperiodic system. Without 
disorder (solid and solid dotted lines) and with disorder (dashed and dashed dotted 
lines). For j3 > 0.5 the response j(t) is positive but for (5 < 0.5 the response is 
positive at small times and negative at long times (backscattering) . 

Metal-Insulator transition 

Let us discuss now the role of quantum interferences according to the scaling 
theory of localization. As explained in part 3.1 the central quantity is the 
conductance of a cube with a size equal to the elastic mean free path L(E F , r). 



g ~ e 2 n(E F )D(E F , r)L(E F , r) (139) 



Since the quantum interferences effect have not the possibility to operate at 
smaller length scale than L(E F , r) then this quantity can be computed with 
the RTA. 

We assume that AXq (E F ,t) is nearly equal to the asymptotic form At 213 then 
the typical propagation length L(E F , r) on a time scale r is of the order of: 



L 2 (E F , t) ~ At 213 (140) 
and the diffusivity is after (105) of the order of: 



D{E f ,t)~^t 2 P- 1 (141) 

From this one deduces that the conductance of a cube of size L(E F , r) is given 
approximately by 



g ~ e 2^£l A 3/2 r 3/3-l (M2) 
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From this expression one concludes that if (5 < 1/3 the conductance tends 
to zero at large r. This means that the system becomes insulating when the 
disorder decreases. This is what happens in the case of crystals if the Non 
Boltzmann contribution to transport dominates (see part 3.3). 

Yet one must note that due to the Guarneri inequality [6] the spectrum be- 
comes singular continuous for a three dimensional system with f3 < 1/3. In 
that case the density of states cannot be considered as a constant in the per- 
fect system. Thus we have to assume that the density of states is sufficiently 
smooth when averaged on the energy scale given by the inverse scattering time 
1/r. 



4 Evidence of anomalous quantum diffusion in quasicrystals and 
approximants 

In this part we present briefly the experimental transport properties of phases 
such as AlMnSi, AlPdMn and AlCuFe or AlPdRe. These experimental trans- 
port properties indicate a conduction mode which is neither metallic nor semi- 
conducting. For the a-AlMnSi phase, recent ab-initio computations are pre- 
sented, which confirm the existence of an anomalous diffusion and allow for a 
semi-quantitative ab-initio computation of conductivity. Concerning AlCuFe 
and related quasiperiodic phases, which cannot be addressed by band struc- 
ture calculations, we present a phenomenological model. This model based on 
anomalous quantum diffusion provides a coherent interpretation of the strange 
electronic transport of these systems. 

4- 1 Experimental transport properties of icosahedral and related approximant 
phases 

Quasicrystals of high structural quality reveal unusual transport properties 
[21,22,23,25] (figure 16). For instance, one of the main features is the low 
conductivity a iK = 100 — 200 ficm" 1 for icosahedral AlPdMn and AlCuFe 
and a 4K < 1 Ocm' 1 for AlPdRe [2,3,4,26,27,28,29], although the DOS still 
has a metallic character. This means that the high resistivity is due mainly to 
a small diffusivity of electrons. 

Experimental measurements show that approximants phases like a-AlMnSi [23], 
1/1 AlCuFeSi [30], R-AlCuFe [23], 1/1 AlReSi [31,32] etc., have transport 
properties similar to those of quasicrystals AlPdMn and AlCuFe. This sug- 
gests that the local atomic order on the length scale of the unit cell, i.e. 
10 — 30 A, determines their transport properties. As atomic medium-range or- 
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Fig. 16. Order of magnitude and schematic temperature dependencies of the re- 
sistivity of icosahedral quasicrystals compared to amorphous and metallic crystals. 
From [25]. 

der of quasicrystals and approximants are similar, it should also be important 
in the understanding of transport properties of quasicrystals. This remark is 
confirmed by the fact that A1TM crystals with a small unit cell (typically 
less than 50 atoms in a unit cell) does not exhibit such particular transport 
properties. 

The resistivity, p = 1/er, of crystals with a small unit cell increases as temper- 
ature T increases and generally follows the Mathiessen rule: 



p(T)= Po + Ap(T), 



(143) 



By contrast, the resistivity of some quasicrystals and approximants (AlPdMn, 
AlCuFe) decreases when temperature increases , and their conductivity follows 
approximatively an "inverse Mathiessen rule" [33,23]: 

a(T) = (To + Ao-(T). (144) 



Besides, after annealing sample, with a consequent reduction of the structural 
defects, the resistivity of quasicrystals and approximants increases. The re- 
lation between the particular transport properties of these phases and their 
structure is still debated. For AlPdMn quasicrystals, J.J. Prejean et al. [34] 
found that local defects might be related with the occurrence of Mn atoms with 
localized magnetic moment. Thus, magnetic properties, transport properties 
and structural quality are intimately linked for those complex phases. 



44 



Another remarkable experimental result is the linear energy dependence of the 
optical conductivity of AlCuFe and the absence of Drude peak [35,36]. 

The icosahedral AlPdRe is the most resistive known quasicrystalline material 
[27,28]. This material displays a strong decrease of the conductivity when the 
temperature is reduced and the conductivity value can fall below 1 (ficm) -1 
at 4K. Although the behavior depends strongly on the composition and the 
preparation of the sample, many authors [2,3,4,26,27,28,29] have reported that 
AlPdRe quasicrystal are very close to the metal-insulator transition. Three 
successive regimes are revealed [28] as the temperature is increased to room 
temperature: a low temperature variable range hopping-like behavior, followed 
by a Thouless regime and a high temperature critical regime. 

Experimentally a low density of states (DOS) at the Fermi energy Ep is usually 
measured in quasicrystals and their crystalline approximants. For instance, a 
density of states at E-p reduced by ~ 1/3 from its free electrons value is 
measured in i-AlCuLi and R-AlLiCu approximant [22,37,38,39]. The presence 
of the pseudogap in these phases is confirmed by photo-emission measure- 
ments [40] and NMR experiments [41]. 

For icosahedral phases containing transition metal (TM) elements, specific 
heat measurement indicate a DOS at E-p of ~ 1/3 of the free electron value 
for i-AlCuFe and ~ 1/10 for i-AlCuRu [42] and i-AlPdRe [43]. From photo- 
emission spectroscopy the pseudogap in the DOS is confirmed for many icosa- 
hedral quasicrystals in the systems: AIMn (metastable) [44], AlMnSi [44,45], 
AlCuFe [46,47,48,49,50,51], AlCuFeCr [49], AlPdMn [52,53,54,55], AlCuRu [56], 
AlPdRe [54] . The pseudogap has been also measured in many approximants of 
quasicrystals. For instance R- AlCuFe [41,23], 1/1 AlCuFeSi [30] a-AlMnSi [23,45], 
1/1 AlCuRuSi [57,58], 1/1 AlReSi [32] have a DOS at E-p reduced by a similar 
factor as in i-AlCuTM and i-AlPdMn. A pseudogap near Ep has been also 
confirmed by ab-initio calculations of the electronic structure in many icosa- 
hedral approximants (figure 17, see below and see also the chapter by Ishii 
and Fujiwara in this book). 

It has also been shown experimentally that transition metal elements have a 
important role on the unusual transport properties of quasicrystals and related 
phases [59,60,61,62,63,64]. 

4-2 Ab-initio electronic structure and quantum diffusion in perfect approxi- 
mants 

Density of states 

Electronic structure determinations have been performed in the frame-work of 
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Fig. 17. LMTO total DOS of Al 6 Mn [63] and a-Al 69 . 6 Sii3.oMni 7 .4. 

density functional formalism in the local density approximation (LDA) by us- 
ing the self-consistent Tight-Binding (TB) Linear Muffin Tin Orbital (LMTO) 
method in the Atomic Sphere Approximation (ASA) [65]. 

The LMTO DOS of an a-AIMn idealized approximant (Elser-Henley model 
[66,67]) has been first calculated by T. Fujiwara [69,70]. This original work 
shows the presence of a Hume-Rothery pseudo-gap near the Fermi energy, 
Ep, in agreement with experimental results [22,23] (see also figure 17). The 
role of the transition metal (TM) element in the pseudo-gap formation has 
been shown from ab-initio calculations [74] and experiments [63]. Indeed the 
formation of the pseudo-gap results from a strong sp-d coupling associated 
to an ordered sub-lattice of TM atoms. Just as for Hume-Rothery phases a 
description of the band energy can be made in terms of pair interactions. 
It has been shown that a medium-range Mn-Mn interaction mediated by 
the sp(Al)-d(Mn) hybridization plays a determinant role in the occurrence 
of the pseudo-gap [68,69,70,71,72,73,74,75,76,77]. It is thus essential to take 
into account the chemical nature of the elements to analyze the electronic 
properties of approximants. The electronic structures of simpler crystals such 
as orthorhombic Al 6 Mn, cubic Al 12 Mn, present [74] also a pseudo-gap near 
E F which is less pronounced than in complex approximants phases. E.S. Zijl- 
stra and S.K. Bose [80] show that Si atoms are in substitution with some Al 
atoms in the a-phase. The main effect of Si is to shift E F in the pseudo-gap 
in agreement with the Hume-Rothery mechanism that minimizes the band 
energy. 
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Fig. 18. Variation An(E) of the DOS due to 1 Mn atom (Virtual Bound State), 1 
Mn icosahedron, and 1 icosahedron of Mn icosahedra obtained after an inflation by 
a factor (p 2 of an initial Mn icosahedron. E-p ~ Ed = 10.88 eV. From [85]. 

Role of clusters 

As for the local atomic order, one of the characteristics of the quasicrystals 
and approximants, is the occurrence of atomic clusters on a scale of 10-30 
A [81]. Nevertheless the clusters are not well defined because some of them 
overlap, and in addition there are a lot of glue atoms. The role of clusters 
has been much debated in particular by C. Janot [84] and G. Trambly de 
Laissardiere [85]. C. Janot considers as a reference clusters that are isolated 
in vacuum but it is more realistic to consider a model of clusters that are not 
isolated but are embedded in metallic medium. The model [85,86] is based on a 
standard description of intermetallic alloys. Considering the cluster embedded 
in a metallic medium, the variation An(E) of the DOS due to the cluster is 
calculated. For electrons, which have energy in the vicinity of the Fermi level, 
transition atoms (such as Mn and Fe) are strong scatters whereas Al atoms 
are weak scatters. Then, following a standard approximation, the potential of 
Al atoms was neglected in reference [85]. 

In the figure 18, An(E) due to different clusters are shown. The Mn icosa- 
hedron is the actual Mn icosahedron of the a-AlMnSi approximant. As an 
example of a larger cluster, we consider one icosahedron of Mn icosahedra, 
which appeared in the structural model proposed by C. Janot [84]. 

An(E) of clusters exhibits strong deviations from the Virtual Bound States 
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(1 Mn atom) [87,88]. Indeed several peaks and shoulders appear. The width 
of the most narrow peaks (50 — 100 meV) are comparable to the fine peaks of 
the calculated DOS in the approximants. Each peak indicates a resonance due 
to the scattering by the cluster. These peaks correspond to states "localized" 
by the icosahedron or the icosahedron of icosahedra. They are not eigenstate, 
they have finite lifetime of the order of h/SE, where SE is the width of the 
peak. Therefore, the stronger the effect of the localization by cluster is, the 
narrower is the peak. A large lifetime is the proof of a localization, but in 
the real space these states have a quite large extension on length scale of the 
cluster or the cluster of clusters. 

The physical origin of these states can be understood as follows. Let us con- 
sider incident electrons, with energy E closed to E-p, scattered by the cluster. 
In AIMn alloys E F ~ E d) where E d is the energy of the d-orbital. In this 
energy range, the potential of the Mn atom is strong and the Mn atoms can 
roughly be considered as hard spheres with radius of the order of the d-orbital 
size (~ 0.5 A). By an effect similar to that of a Faraday cage, electrons can 
by confined by the cluster provided that their wavelength A satisfies A > /, 
where / is the distance between two hard spheres. In the case of a-AlMnSi 
approximant, A ~ 3.6 A (if we assume a free electron band and E F = 10.33 
eV) and the distance / is about 3.8 A. Consequently, we expect to observe such 
a confinement. This effect is a multiple scattering effect, and it is not due to 
an overlap between d-orbitals because Mn atoms are not first neighbor. We 
have also shown that these resonances are very sensitive to the geometry of 
the cluster [86]. For instance, they disappear quickly when the radius of the 
Mn icosahedron increases. 

Quantum diffusion in perfect crystals 

In the following we present calculation of the quantum diffusion in perfect 
crystalline systems. Some works have already been done from ab-initio cal- 
culation and give indication of none ballistic diffusion [17,89,90]. We consider 
the a-AlMnSi approximant and compare it with simpler crystals orthorhom- 
bic Al 6 Mn, and cubic Al^Mn [91,92,93]. For the a-AlMnSi phase, we use the 
experimental atomic structure [94] and the Si positions proposed by Ref. [80] 
with composition: a-Al69.6Sii3.oMni7.4. In figure 19, the total DOS n of the 
a-AlMnSi phase is presented versus the energy. The total density of states is 
characterized [69] by a pseudogap near the Fermi energy E-p. Following the 
Hume-Rothery condition, it is expected that the most realistic value of E F 
in the actual a-phase corresponds to the minimum of the pseudo-gap, i.e. 
E-p — i?F(LMTO) = — 0.163 eV for our calculation. 

We compute the velocity correlation function C(E,t) for crystals (complex 
approximants and simple crystals). In equations (12), (37), the average is 

— * — * 

obtained by taking the eigenstates for each k vector with and energy E n {k) 
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Fig. 19. (A) LMTO total DOS n, (B) conductivity a, and (C) diffusivity D, in the 
cubic approximant a-Al69.6Sii3.oMni7.4. From [91]. 

such as 



E - -AE < E n (k) <E + -AE. 
2 2 



(145) 



AE is the energy resolution of the calculation. The calculated C(E,t) is sen- 

— * 

sitive to the number N k of k vectors in the first Brillouin zone when N k is too 
small. Therefore N k is increased until C(E,t) does not depend significantly 
on N k . For Al, Ali 2 Mn and a-Al 69 . 6 Sii3.oMni 7 . 4 , AE is equal to 0.272, 0.272, 
and 0.0272 eV, respectively, and N k is equal to 120 3 , 40 3 and 32 3 , respectively. 
C(E F ,t) for the cubic approximant a-AlMnSi is shown in figure 20. 

After (12) and (37), C(E F , t) is the sum of a constant Boltzmann term Cb(E, t) 
and a non Boltzmann term containing oscillating terms that average to zero 
on long time scale: 



C(E,t) = C B (E,t) + C NB (E,t) 



(146) 
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Fig. 20. Velocity correlation function C(E F ,t) (m 2 s~ 2 ) in a-Al69.6Sii3.oMni7.4 ver- 
sus large time t The dashed lines are the corresponding Boltzmann velocity corre- 
lation function Cb(-Ef,0 = 2-uf. From [92]. 



C B (E,t)=2{V 3 



(147) 



Urn j C NB (E, t) e~*/ r dt = (148) 
o 

where is the square of Boltzmann velocity (intra-band velocity) along the X 
direction at the Fermi energy: v F = 9.4 10 7 , 3.6 10 7 , and 2.7 10 6 cm.s" 1 , for Al, 
Ali 2 Mn and a-Al 69 . 6 Sii3.oMni7.4, respectively. This last result is very similar 
to the original work of T. Fujiwara et al. for the a-Alii 4 Mn 2 4 (with the atomic 
structure model of Elser- Henley) [70], for a model of icosahedral approximant 
AlCuFe [82]. The strong reduction of v F in the approximant phase with respect 
to simple crystal phases shows the importance of a quasiperiodic medium- 
range order (up to distances equal to 12-20 A). This leads to a very small 
Boltzmann conductivity for approximants [70,82]. In the case of a decagonal 
approximant AlCuCo, a strong anisotropy has been found between v F in the 
"pseudo" quasiperiodic directions and vp in the periodic direction [83]. 

On small time scale t (figure 21), C(E F ,t) and CB(E F ,t) differ, and there is 
a new difference between approximant and simple crystal. In the case of Al 
(f.c.c.) phase, C(E F ,t) is always positive, and the Boltzmann value is reached 
rapidly when t increases. But for some t values the velocity correlation function 
C(E F ,t) is negative for Al 12 Mn and a-Al 114 Mn 24 . That means that at these 
times the phenomenon of backs cattering occurs. 

The transport properties depend on the average value of C(E F ,t) on a time 
scale equals to the scattering time r [6,97] (see for instance equation (29)). A 
realistic value of r has been estimated to about 10 -14 s [33]. For the simple crys- 
tals Al 12 Mn, C(E F ,t) is mainly positive when t > 2 10~ 15 s. But for the com- 
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Fig. 21. Velocity correlation function C(Ep,t) (m 2 s -2 ) versus small time t 
The dashed lines are the corresponding Boltzmann velocity correlation function 
C B (E F ,t) = 2v|. From [92]. 

plex approximant a-Al 114 Mn24, a lot of t values correspond to C(Ep,t) < 0, 
even when t is close to r or larger (figure 20). Therefore, in the case of Al 12 Mn, 
the backscattering (negative value of C(Ep,t)) should have a negligible effect 
on the transport properties, whereas this effect must be determinant for the 
approximant. 

As discussed in part 3 the phenomenon of backscattering is associated to 
unusual quantum diffusion. It is illustrated on the plot of the average spreading 
of states AX 2 versus time t (figure 22). After (40), AX 2 results in two term: 

AX 2 (E, t) = V B (E) 2 t 2 + AX 2 B (E, t), (149) 

A Boltzmann term V B (E) 2 t 2 and a non-Boltzmann term. The non-Boltzmann 
contribution, AX^ B , which comes from the non-diagonal matrix element (44), 
increases very rapidly and saturates to a maximum value of the order of the 
square size of the unit cell. In the a- approximant, at small time t, AXg is 
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Fig. 22. Square spreading AX 2 of electrons states with Fermi energy Ep versus 
time t, in the cubic approximant a-Al6g.6Sii3.oMni7.4. AX 2 = AX^ + AX^ B (see 
text). From [91]. 

smaller than in Al due to a very small velocity Vp of the electron with energy 
Ep. The calculated Vp is equal to 2.7 10 7 cm.s~ 1 , which is about 30 times 
smaller than aluminum values. Thus a-AlMnSi is a non-conventional metal at 
these time scale i.e. when the scattering time is t < t* where r* ~ 3 1CT 14 s. In 
a normal crystal, the AX^ B (t) term is negligible with respect to the Boltzmann 
term AX B (t). On the contrary, in the approximant both terms have the same 
magnitude at the realistic scattering times scale, typically a fews lCT 14 s. 

4-3 Ab-initio RTA model for the conductivity of approximants 

Within a relaxation time approximation the diffusivity D(E,u) is calculated. 
At low frequency one gets 



Re a(E, to) = e 2 n(E)D(E, to) 



(150) 



D(E, uj) = D B (E, u) + D NB (E, u) 



(151) 



D B (E,cu) 



V 2 



T 



(152) 



and 




(153) 
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Fig. 23. Ab-initio dc conductivity a in cubic approximant a-Al69.6Sii3.oMni7.4 ver- 
sus inverse scattering time. 

The D B values for a-Al 69 6 Sii3.oMn 17 . 4 (figure 19) are similar in magnitude 
to those obtained by T. Fujiwara et al. [69] for the idealized approximant a- 
Aln4Mn24 approximant. -Dnb is almost independent on E, whereas the 
values depend strongly on E and is particularly small in the pseudo-gap. 

The predicted static conductivity (dc conductivity) of the a-AlMnSi phase, 
assuming the value of the Fermi energy given above i.e. Ep — -E'f(lmto) — 
— 0.163 eV, is shown figures 23 and 24 versus the inverse scattering time. Two 
regimes appear clearly: the metallic regime (Boltzmann regime) at large scat- 
tering time, t > t*, and the insulating like regime (non Boltzmann regime) 
at small scattering time, r < r*. r* = 3.071CT 14 s is defined as the time 
for which the Boltzmann and non-Boltzmann contributions are equal. As ex- 
pected from our model, <tnb is almost proportional to 1/r. Therefore, in the 
non Boltzmann regime, the conductivity increases with disorder as observed 
experimentally. For realistic r values, t < r* [33], onb dominates and a in- 
creases when 1/r increases i.e. when defects or temperature increases, a varies 
from 250 (ft cm)" 1 for r = 3.310~ 14 s, to 2000 (ficm)- 1 for r = 10~ 15 s. This 
is consistent with experimental results in a-AlMnSi: er(4 K) ~ 200 (ficm) -1 
and <r(900 K) ~ 2000 (Ocm) -1 and with standard estimates for the scatter- 
ing time in these systems [23]. Furthermore for r equal to a few 10~ 14 s, i.e. 
when the Boltzmann term is negligible, the mean free path is given by the 
square root of the saturation value of AX^ B and is of the order of 15 A. This 
is in agreement with estimates in the literature [23]. As discussed in part 3 
this means also that the systems is far from the Anderson transition despite 
its low conductivity. From the ab-initio calculations the estimated value of 
the ratio go/g c fc> r the a-AlMnSi phase is about 2 — 3. This means that this 
system is always metallic as discussed in part 3.3. According to figure 13 the 
a-AlMnSi phase is a phase of type (a). 

Optical conductivity 
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Fig. 24. Ab-initio electrical conductivity a(Ep) (logarithmic scale) versus 
inverse scattering time 1/r (logarithmic scale). (A) In cubic approximant 
a-Al69.6Sii3.oMni7.4: a(Ep) is the sum of a ballistic term (Boltzmann term), 
0"b = e n(Ep)VpT, and a non ballistic term (non Boltzmann term), <tnb- (B) In 
pure Al (f.c.c), the Boltzmann term dominates, and the model is compatible with 
a simple Drude model (dashed line). In cubic Al^Mn and orthorhombic A^Mn 
crystal, the model predicts also a Boltzmann behavior as expected experimentally. 
From [91]. 
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Fig. 25. Inter-band optical conductivity, o"nb(^), hi oAlii4Mn24 calculated from 
the LMTO results, with a relaxation time r equals to infinity. From T. Fujiwara et 
al. [70]. 
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Fig. 26. Ab-initio optical conductivity ct(uj) in cubic approximant 
a-Al69.6Sii3.oMni7.4 for three r values, uj is the pulsation. For r = 3r*, the 
non Boltzmann conductivity <tnb is smaller than Lorentzian of the Boltzmann 
conductivity cr B , cr B (w) = cr B (0)/(l + u 2 t 2 ). For r = r* = 3.07 x 10~ 14 s, 
cnb(0) = <T B (0). For r = r*/3, the non Boltzmann conductivity cjn B (w) dominates. 

Within the relaxation time approximation used here, the optical conductivity 
<j{uj) is the sum of two terms. The Boltzmann contribution (<r B (u;), diagonal 
elements of the velocity operator) gives rise to the so-called Drude peak and 
the non Boltzmann conductivity gives rise to a nearly frequency independent 
contribution. This is a consequence of the fact that AX^ B (E,t) is nearly 
constant on the time scale of r (see part 3). 



T. Fujiwara et al. [70] has also estimated the optical conductivity from the 
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Fig. 27. Ab-initio dc conductivity a in an hypothetical cubic approximant 
a-Al69.6Sii3.oCui7.4 versus inverse scattering time. 



LMTO band dispersion of a a-Al 114 Mn 2 4 (figure 25). This calculation repro- 
duces the linearity and the peak position observed experimentally. Our ab- 
initio calculation (figure 26) confirms that a Drude peak can be identified in 
the Boltzmann regime, r > r*, whereas in the non Boltzmann regime, r < r*, 
the Drude peak disappears. 



The role of transition metal elements (TM = Fe, Mn, Co, Pd, Re) in the elec- 
tronic structure of quasicrystals and related phases as been often discussed 
in the literature [68,69,70,71,72,73,74,75,76,77,78,79]. Because of their strong 
scattering potential with respect to Al(Si) atoms, TM elements play a cru- 
cial role in the formation of the Hume-Rothery pseudogap that contributes 
to the stability of these phases. This effect is related to an effective medium 
range interaction between TM atoms mediated by a strong sp(Al)-d(TM) hy- 
bridization [74] . TM elements have also a very important role on the transport 
properties. As an example, it is shown in the previous paragraph how a Mn- 
cluster can "localize" electrons [85,86]. 



To evaluate the effect of TM elements on the conductivity calculation in the 
RTA, we have considered an hypothetical a-Al 69 . 6 Sii3.oCui 7 .4 constructed by 
putting Cu atoms in place of Mn atoms in the actual a-Al 69 . 6 Sii3.oMn 17 . 4 
structure. Cu atoms have almost the same number of sp electrons as Mn 
atoms, but their d DOS is very small at Ep. Therefore in a-Al 69 6 Sii3.oCui7.4, 
the effect of sp(Al)-d(TM) hybridization on electronic states with energy near 
Ep is very small. As a result, the pseudogap disappears in total DOS, and the 
dc-conductivity is now ballistic (metallic) as shown on figure 27. 
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4-4 Phenomenological model for the low frequency conductivity of AlCuFe 
quasicrystals 

The ab-initio calculations which rest on the Bloch theorem are applicable to 
approximants only. Here we present a phenomenological model of the optical 
conductivity of AlCuFe QCs, which should be approximately valid also for 
the related QC phases AlPdMn or AlFeCr [95,96]. This anomalous diffusion 
model allows to derive an analytical expression for the conductivity that fits 
experiments very well. In particular the model explains quantitatively the 
main experimental facts: 

- the increase of conductivity with disorder 

- the "inverse Mathiessen rule" [23,34] that is the fact that the increases of 
conductivity due to different sources of scattering are additive 

- the absence of the Drude peak 

One needs first a model of conductivity of the perfect system at all frequencies. 
For low frequencies, according to the discussion in part 2 we assume: 

/ 1 i \ 1_2 ' a 



A small value of /3 {(3 <C 1 ) is imposed by the nearly linear experimental vari- 
ation of Re <t(u) at uo < 8000 cur 1 (see figure 28). This means that the system 
without defects would be insulating. In (154) we take f3 — 0, u\ ~ 8000 cm" 1 
and o" ~ 6000 (f2cm) -1 . At higher frequencies one uses other analytical ex- 
pressions. For 8000 cm -1 < uo < 25 000 cm -1 a polynomial of uo reproduces 
the experimental value. For uo > 25 000 cm -1 , we take the Drude expression 
according to [35]. Let us note that the experimental uncertainty on the high 
frequency conductivity [35] has essentially no effect on the results presented 
here. 

Within the RTA (3,25) one has for the optical conductivity: 



i.e. the real part Reff(w,r) of the conductivity of the system with defects is 
the convolution of Re a (u) of the perfect system and of a Lorentzian of width 





(154) 




(155) 



1/r. 
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Fig. 28. Real part of the conductivity as a function of the frequency for differ- 
ent cases. Line with circles: experimental conductivity of an AlCuFe QC [9]. Thin 
line: conductivity of the model without defect Re<7o(w) = ai{\u\/ui). Thick line: 
conductivity of the model with defect, for r = 3 10~ 14 s. From [97]. 

As shown in part 2 for uj < cu 1 the conductivity is well represented by 



A 

Re a ~ — 

T 



a + log 



UJ\T 



+ ujt Arctg(cjr) 



(156) 



where A = 2a /nuJiThe analytical expression (94) with a ~ 0.7, fauJi ~ 1 eV 
and (Jx ~ 6000 (11cm)" 1 describes well the electronic conductivity in figures 
28, 29. 

Let us focus on the low frequency conductivity [uj < uji) which is the real test 
of the model (155,154,156). Figure 28 gives a comparison of the experimental 
Rea exp (uj) for AlCuFe [35] with the theoretical Rea(uj,T). The scattering time 
r is chosen to reproduce the experimental dc-conductivity a dc — 350 (ficm) -1 . 
One finds r ~ 310 _14 s which is rather long, in agreement with the high 
structural quality of these systems. The fit is good except for the peak in the 
experimental curve around 200 cm -1 . This peak is attributed to the conduc- 
tivity of phonons [35] which is not incorporated in the model. The mean-free 
path A is related to the scattering time r and to the diffusivity D through 
D = A 2 /3t. One estimates [23,34] D > 0.2cm 2 /s and since r ~ 3 10~ 14 s, one 
has A > 15 — 20 A. From (155) and (156) one gets the dc-conductivity tTd c (r) 
as a function of the relaxation time r (see figure 29). o" rfc (r) increases with 
1/r and varies nearly linearly with 1/r on a large range of values of o- dc {r) 
i.e. a dc ~ A + B/t. For two independent sources of scattering characterized 
by scattering times T\ and T2 it is common that the inverse relaxation times 
add. Then 1/r ~ 1/n + l/r 2 and a dc ~ A + B/t ~ A + B/n + B/t 2 . Thus 
each source of disorder gives its contribution to the conductivity in agreement 
with the "inverse Mathiessen rule" [23]. 
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Fig. 29. Thick line: Variation of Od c with x = 1000/cjit. t is given by 
r = (6.6/x) 10~ 13 s. The straight thin line shows that a^e varies nearly linearly 
with 1/t in the range a dc = 150 700 (ficm) -1 . From [97]. 
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Fig. 30. Value of j(t) deduced from the experimental conductivity. The negative 
value of j(t) at large times indicates backscattering. From [97]. 

The present phenomenological model treats the disorder within the relaxation 
time approximation (RTA). Indeed, as shown now, the RTA is applicable to 
AlCuFe QCs, at least for T < 200 — 300 K. A first indication is that quantum 
interferences have been found for T < 200 — 300 K [23] . They indicate that the 
main scattering sources are elastic in this temperature range. Indeed, if the 
dominant scattering were inelastic the coherence of the electron wavefunction 
would be lost at each scattering event. In that case there would be no inter- 
ferences in the diffusive regime. In addition the experimental fits [23,34] show 
that the quantum interferences and the electron-electron interaction give only 
a correction to the conductivity. Therefore, as the elastic scattering dominates 
and as quantum interferences are weak, the RTA is a good approximation for 
the AlCuFe QC studied in [35] at least at T < 200 - 300^. In particular 
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a scenario of hopping between localized critical states, such as proposed by 
Janot [84] is not consistent with the present analysis. 

Note that the model is consistent with the observed weak-localization effects. 
Indeed in the context of the scaling theory of localization [98] the importance 
of quantum interferences depends on the ratio between the dc- conductivity 
of the system a dc and the Mott value OMott — 600 (f2cm) _1 /A where A is the 
mean-free path expressed in Angstroms. If R = cr dc /cr Mott ^> 1 the effect of 
the quantum interferences on a dc is small. Here R = cr dc / aMott — 5 — 10 and 
the localization effects are only corrections. 

To conclude, the dynamics of electrons in AlCuFe quaiscrystals and related 
systems such as AlPdMn [95], AlCrFe [96] quasicrystals is not free electron 
like. The minimum of optical conductivity at low frequency and the increase 
of dc-conductivity with disorder (in the RTA scheme) are intimately related 
to the backscattering (figure 30) or equivalently to an anomalous diffusion 
law with (3 < 0.5 (Ref. [6]). The low value of j3 implies a slow anomalous 
diffusion L(t) oc t 13 , which suggests the proximity to a localized state. This 
value P ~ is consistent with the results on approximant phases. Indeed the 
Non Boltzmann term in the quantum diffusion AX^ B (t) which saturates very 
quickly is analogous to a quantum diffusion with (3 ~ 0. 



5 Conclusion 

In this chapter we concentrated on quantum diffusion and electronic conduc- 
tion properties in quasiperiodic and periodic systems. We found that devi- 
ations from the standard ballistic propagation exist either in quasiperiodic 
or even in periodic systems. This anomalous diffusion mode has deep conse- 
quences on the conduction properties at zero and low frequency. 

The anomalous diffusion mode is related to a tendency to localization and to a 
phenomenon of backscattering which is well known in disordered systems. The 
phenomenon of backscattering is the fact that an impulse of electric field cre- 
ates a current density J(t) which is opposite to the electric field at large time. 
Backscattering is associated with an increase of conductivity with frequency 
and disorder. 

The physics of phonons in quasicrystals could be affected by the anomalous 
diffusion phenomenon. In particular it has been argued that the heat conduc- 
tivity could be sensitive to this effect [99] . 

The concepts developed here open also a new insight in the physics of cor- 
related systems. Indeed recent studies of some heavy fermions or polaronic 
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systems [100,101,102], where charge carriers are also slow, show that their 
conduction properties present a deep analogy with those described here. In 
particular a transition from a metallic like regime at low temperature where 
scattering is weak to an insulating like regime at higher temperature with a 
stronger scattering is observed. 
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